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ABSTRACT 

!c t* ibo*rt :tut a wide tUm of setcMial problem* involving aniso::op?c media eaa be ;ra.u» 
forced into «?Ti:v*jn;{ problems Evolving oaiy isotropic tr.cdi*. By means .jf svica *rar.iforn;irir,ii 

i< is pOSWWc, tft A U/M numb<f Of <we». to determine! tbe apparent c**i*li vitle» *bic!> iv-intf be 

observed la anisotropic ionwtiOWk urn* riectr<*le-:yp« resistivity Sow** devices. Di*x»inr, is 
tfvca of *n ilOflite, anisotropic medium *i;.S» *rxl without borehate, ft! two f*mi-in4flite j.iij«tra?ie 
b«ds fwitbejut b«r*a$le]l. »0d of l rtun wiropie bed bouaeeri by anj^tr^pic u)!iec:it :crjr.iiir,^s 
»«*itaout Sesehoiei. An-irUerprtU'-ioci cb*r; tar :he nor:naJ dtvk= la proer.ted for thick. :icn>tn vac' n\. 
*ai»t:optc bees pcaetraied by a wwtaolc. 



INTRODUCTION 

In many sedimentary strata, electric current flows more easily in a direction 
parallel to the bedding planes than transversely lo them. The reason ioi this 
anisotropy is chat a great number of mineral crystals possess a flat or elongated 
shape (e.jr., mica and kaolin!. At the time they were laid down, they natuialiy 
took an orientation parallel to the plane of sedimentation. The Interstices in ;he 
formations are, therefore, parallel to the bedding plane, and the current is able 
to. travel with greater facility along these interstices which generally contain 
mineralized water. Such electrical anisotropy is observed mostly in shales. 

If a cylindrical sampic is cut from a formation, parallel to the bedding p/anes, 
the resistivity of this sample measured with a current Sowing along its axis is 
called the longitudinal resistivity R«- If a similar cylinder is cm perpendicular ty 
the bedding, the resistivity measured with a current flowing along its axis is called 
the t ransver sal resistivity Rv. The anisotropy coefficient X, by definition, is cqui? 
to ^Rv/Rh- Laboratory measurement* have shown that \ may range from I :o 
about 2.5 In different shales. 

Furthermore, tbe formations *re often made up of a series of relatively thin 
beds having different lithologic diaracteristics and, therefore, different resitiv- 
itiea (as, for example, sequences of thin shales and hard streaks). In electrode 
system* 'used tor electrical prospecting and ior electrical logging, the distances 
between the electrodes are great enough that the volume involved in a areas;,--." 
rnent may include several such t^in beds. Since* in this situation, the cur:*"' 
3ow» more easily along the more conductive streaks than transversely t:> 
series o£ beds, there i* effective anisotropy. The effects on resistivity mea»u.'L- 
ments of this ■■macroscopic" anisotropy an? added to the effects of the ar.iso'r..; y 
due to the microscopic structure of the sediments which was described above 

The cicctrtcal anisotropy of sedimentary formations was recognised *> 

■ M»pu*sipt received by the EdHar January 18. 

• Seblumbergti W*U Surveying Corporation. Rj^ictDcUi, Ccaa=eiio«. 
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FOXMATIQP ASISOT&OPY EFFECTS OX XEStSTlVITY MEASUREMENTS 

C. Schlumberger ;l°2Q> early in the history of electrical prospecting. Later, in ■ 
1932, a method of computing the electric field in anisotropic media was presented 
by K. Maiilet and H. G. Doll. Practical use Jias aiso been made of the aniso- 
tropic properites of sediments for the determination of the directions of the 
strike and the dip. both in surface prospecting and in well logging (Schlumberger 
and Scblumbcrscr, 1950; Schlumberger, Schlumberger, and Leonardos, 
Schlumbergcr, Scblumbetger r and Doll, 1923). 

One important result brought to light in the papers by Schlumberger. Schlum- 
berger, and Leonardon (1934) and by Maiilet and Doll (1932) « that the resistiv- 
ity measured with a system of electrodes aligned perpendicularly ;o the bedding 
planes is equal to the longitudinal resistivity (paradox oi anisotropy). In the par- 
ticular case of ekctrieal logging in boreholes, the measuring device is generally 
vertical, and, whenever the formation dip is small or nil, the measured resistivity 
accordingly is equal to the, longitudinal resistivity, provided the electrode spacing 
is large enough for the effect of the borehole to be negligible. 

Recently, some features observed on resistivity logs have called attention to 
the need for a more complete investigation of the effect of anisotropy. One of these 
phenomena is the difference in readings between the Induction Log 1 and the 14- 
inch Normal 1 at the level of certain shales (the former being systematically 
lover than the latter). Another phenomenon is the abnormally low values read 
v^tb the lateral device 1 at the level of certain sands bounded by shales. 

A mathematical study has accordingly been made of the effect of anisotropy 
on resistivity logs made with electrode devices, with special attention to the 
following cases: 

a) Homogeneous, infinitely thick, anisotropic formations, Ix&verscd by a 
borehole; 

b) Sequences of horizontal* homogeneous, anisotropic beds, borehole tffect 
neglected. 

It will be shown that the evaluation of the anisotropy effect is possible by sub- 
stituting for each anisotropic formation an appropriate isotropic medium. Using 
such equivalences, the responses of normal and lateral devices in anisotropic 
media can often be determined, 

Departure curves, giving the apparent resistivity read by a normal device in 
a thick anisotropic bed as a function of the electrode spacing, are presented 
ior several values of the anisotropy coefficient. These curves take into account 

i explained in detail in D«C (1049), in measuneowt* of wautivity by Induction Lowing, 
current loop* eneifdin 5 the &xi* of the bole are induced in the earth, 5iace, for » vertical driJl hole, in 
the ihsencc of formation dip, these currrats <5rculate horiaontaiJy, their flow » affected. only by the 
longitudinal reaiatmty; ana the reading so found » a <bw« of Ka- 

* For Oie *rr»ngcn»enL of the wrcnt and metering «Wcm*Je« iti normal and lateral reriwivrcy 
logging devices, vx Fifiur t !. For a ax« detailed txplatu-t^, jee Chapter I. Section?, ct the Sefalna- 
bergrr Corponrioa'* UurprcMim 8o*f-$o*k i<* Stiiiiixiiy L»tt. A Icj-inch formal L> a normal 

device with rvn AM ifuoof oi 16 uiM- . . f 
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the presence o! a borehole and have been numerically computed using an IBM 
701 Data Processing Machine. 

The behavior of the normal and lateral curves is discussed for the case of i*o 
semi-infinite media with a common boundary, and for the case of a thin. bed. 

OX^EAAL SQVaTXO^S tOK CVttK QX1CAI SYMfCTXY 

On the assumption that all ihe beds are horizontal and the borehole is vertical, 
one may speak of a horizontal ^longitudinal) resistivity R n parallel to the bedding 



SPACING AO 



SPACING AM 

Jl *m 



MIDPOINT 0 



(0) NORMAL (b> LATERAL 

Tic. t. Electrode designations for the limpU&ed norwal and later*] devices referred to in tbi* 
piper, A designates a current electrode, where** ir* sad A* design au« aic**u« electrodes. A current 
rct«f« d«trode, 5, is located *t a tj6« distance Lrwta U» boia device proper. In the idealised Uteri). 
4iaumed in tic copulations. Measure electrode* if and If ire very close together. 

To timpfify the discussion, the electrodes are treated as points. No consideration U mad* of the 
size oC th« electrode* or the presence rf tie insulating mandrel oo wfckb ti*y are mounted. Ti.fi 
eEects ot ticse. aowever, are minor. 



pUaes and of a vertical (transversal) resistivity Rv perpendicular to these planes. 
If the current sources lie on the axis of the borehole, the potential surfaces and 
the line* of current fiow posse** a cylindrical symmetry. 

In a. cylindrical co-ordinate system {». s, $) t Ohm's law may be generate 
to state that the current density J is given in terms of the potential V by the 
equation 

3-3 
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FORMATION AS1S0TR0FV EFFECTS O.V RESIST iYtTY MEASUREMENTS 

where n is a uuit vector parallel to the axis of the borehole (s-axis) and ^ a unit 
vector perpendicular to 2\ and directed away from the borehole axis. 

The nora&ls to the equipotential suriaces, P(p, s)* constant, are parallel to 
the gradient of V, given by 

„ 3K a 7 

* — »i + — *K. (2) 

CO dZ 

Thus, in general, the current flow in an anisotropic formation (Rv^Rft) will not 
be norooal to the potential surfaces. 

The absence of current sources or sinks, everywhere except on the borehole 
axis, is expressed by requiring that the divergence of the current density vanishes 
for »^Q, This leads to the differential equation 

p dp\R„ bp } 6z\Rv *s / w 

For a homogeneous. Isotropic medium, in which «= Rr~ constant, this equation 
reduces to the Laplace equation in cylindrical co-ordinates. 

It is convenient to introduce che m<c»M resistivity R. and the cot£Uit*U of 
aniiotrdpy X, defined by setting 



R =* \/Ra • Rr 



Then 



R u -—> Rr - \R. (5) 

Tilt is assumed that the variation ot the coenldeat oi anlsotropy with p and * 
is expressible in the form 

■*(*), W 

then, by (3) 

By introducing the new variables 

p' = f *nd »' = f W 

^ a /0>> •» • 

and letting 

V(p, s) - V'V, O. W 
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one cads that 

dV 1 dV 
dp i dp' 

and 



<'t0) 



di a ds' 
With these transformations, (.') becomes 

d\"\ 6 / * dV- 



5 — r • Hi- 



J-(JL S 11) + ±(JL?1) = Q . (l2) 

$p \R dp' / dz 4 \R H' / 
Finally, by Introducing a. new value of resistivity at each point, given by 

R* « L R, >\l\ 

p 

one may replace (12) by 

A comparison of equation (14) *'ith (J) shews that in terms of the printed 
co-ordinates (a*, z') the differential equation tor the potential V i% that required 
for an isotropic medium of resistivity R'. Therefore, within the conditions :•«• 
posed by any medium involving anisotropic formations may be replaced by an 
equivalent isotropic one. This is done by mapping each point (p, z) in the ori«:n,t! 
»pace into * corresponding point {»' . of a new space. In this mapping :Sk 
potential remain* the same at corresponding points (tee equation (,&)}, but -.he 
roean resistivity changes from R to R' (-mp'R/fi). U remain* necessary only -.o 
investigate how the boundary conditions are to be transformed. 

TRASSFOftMJLTtON Ot BOU>TDA*Y COWDtTtONS 

The boundary conditions of importance bete are those involving the con- 
tinuity of potential and ot the normal component of current across various 
boundaries. The continuity of potential in the primed system is insured immedi- 
ately by (9), and thus wc need only to consider the transformation of the cur:er.> 
across equivalent sections of a given boundary. 

The current flowing out through the sides of any cylinder of radius ,.whr.« 
a*is is that of the borehole, is given by 

dV 



r * \ dV , 



3 5 
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FOtOiATlOX ASISOTROPY EFFECTS OS R&S1STIVTTY W EASVREUEyTS 

where r, and r t are the : co-ordinates of top and bottom of the cylinder. The cur« •— S 
cent /' flowing out of the corresponding cylinder of radius in the transformed 
spj.es is given by 

r §t ' i w 

1h~-1*p'\ — — di'. (15b) 

J „• J? dp 

Direct substitution from equations t'5). I6i, ',8?, (10), and (13) proves that 

The current flowing in the positive s-dircction through any plane horizontal 
annular surface, bounded by circles of radii pi and Pt, is given by 

which again can be shown to be equal to the current Vy through the correspond- 
ing annular suriace in the primed system. Liing these results, it is easy to show 
;hat the transformation described by equations (8), (9), and leaves un- 
changed the current flowing through corresponding surfaces. In particular, this 
means that the strength of any current source or sink will nor he changed. 

AS INFINITE, HOMOGENEOUS, ANISOTROPIC UEStTC 

For such a medium, X = constant, and one can, by (6), mate one of the follow- '^j 
ing choices; 

r*0») - !, ,(s) - \ (17a) 

or 

f(f) =- X, - I. (17b) 

By equations (8) and (13) the first choice leads to an. expansion in the 
direction given by z'**\z and a resistivity /?' of the equivalent medium equal to 
the mean resistivity R; whereas the second choice corresponds to a radial com- 
pression given by p' = p/X and an equivalent homogeneous resistivity giv«n by . 

The transformation to an isotropic medium represented by the second choice, 
since it involves no change in the electrode spacings, shows at once that theoret- 
ically , if ike borewlc could be nttftti-cd, a normal or lateral device would measure 
Rh, the horizontal resistivity (paradoa of anisotropy). 

For the first choice, which corresponds to a transformation to an isotropic 
medium having a resistivity equal to [he mean resistivity R, it is necessary only 
to take into account the increase by a factor X in the equivalent electrode spac- 
ing. This increase in the spacing will mean that the potential measured by a 
normal sonde, and hence the apparent resistivity measured by it, will be reduced 
by the factor i/7\, and thus one again measures K B — Rf\. For the lawral sonde . 



>e 
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K. S. KUSZ AND J B. XORAH 

the increase in the spacing .40 reduced the signal by 1/X', but the increase in the 
measure electrode span MN iitcxeascs the potential difference between them by 
X, and the boruontal resistivity R* is again measured. 

EFFECT or THE SORXHOtX 
Subsurface resistivity measurements- with an electrode device are generally 
made by towering the device into a drilled bole filled with drilling mud, Trie pres- 
ence oi the electrically conductive fand isotropic) column of drilling mud con- 
taining the device can be taken into account in equation i6j by noting that 

tX* tor p > a, 

where a is the radius of the borebo.le. Also, since X La a function oi p, one sets 
f\0) m X and tii} ■ I; thus, from (8), 

j' - s (19a) 

and 

t p tor P < a \ 

U for P >aJ 
The equivalent isotropic resistivity after the transformation is, from equa- 
tion (13). 

f R m for p < a 

— — — tor p > «, 

U«p' - (Xo - 1)^ 

*here R~i* tbe r«&btivity of the mud filling the hole, and R is the mean resistivity 
of the anisotropic formation. The radial resistivity profile as given by (20) is 
shown in Figure 2. Note that, as the resistivity R' of the equivalent 

isotropic formation approaches asymptotically R» - R/\t. the horizontal resistiv- 
ity. At the edge of the borehole 'J - *) , as seen from (20), R' is equal to J?; but at 

which is less than twice the borehole radius, R'-KR-rRu). At this distance JT 
has returned halfway to its asymptotic value R*. 

From this analysis the departure curves for normal and lateral devices should 
lie above the isotropic curves for a formation of resistivity R B and below those 
for a formation of resistivity A-X^R*. *nd they should approach the former for 
Urge spacing*. This problem is treated quantitatively in the next section, arid.as 
seen in Figure 5, th*** cfc»r»eieri*itc f«aiure* are ind^d borne out. One is \t6, 
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FORMATION AXISOTRQPY EFFECTS OS RESISTIVITY M EASV HEM EST S 

therefore, 10 believe that the transformation to aa equivalent isotropic formation 
»iU ofies serve to gain a qualitative prediction oi the apparent resistivity to be 
expected. 

QCAN'TITATrVE STtTDY Otf TfiS BORtBOtC 

It ha* h«cn shown above that * normal or lateral sonde in an infinite, homo- 
geneous, anisotropic medium will measure Rn and (hai the effect of the borehole 
Through such a medium if to increase the observed apparent resistivity 



R' 



9 

(At 
g 

s 
V 
i 
i 



ft * s/ Ry'R N 
^. ■yR v /ft M 




Fie 2. RatfiaJ plot of the resistivity R' af IB iwtre&te medium equivalent b> a homogencou* 
anisotropic medium pierced by an isotropic boroogpaeeu* uiui tahunn of resistivity fit*. (Plotted for 

as compared to a homogeneous, isotropic formarioa of resistivity In fact* 
the formation will appear as an isotropic medium with a resisdviry lying be 
tween R» and \Ra, depending on Lbe spacing. 

This problem «iU now be formulated quantitatively and some numerical 
results given. Let the. potential in the borehole be V m while that in the formation 
is Vf, Then a solution of (3) roust be found subject to the boundary conditions: 

2R m 

(a) 



(0 



V.,V>-0 



as i — »0, 
IS r— • », 

at p » a, 



(21) 
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xi 2,0 

X'2.5 



Fee }. Normal dtvicc (no lov*»ioa). 
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resulting forinuia has eke same torn as Simpson's formula except that the weights ;• 
:o be assigned :o the ordinate* depend on the magnitude of z*/o- and reduce to 
Simpson's weight* only when ;A/a i* jmali compared to one, h being the width 
of the --intervals in. the numerical integration. 

This method of evaluating the integral in (24)' was programmed for the IBM 
701 Processing Machine and «w found to give satisfactory results even jor 

a highly oscillatory integrand corresponding to z/c » 236. 

The results of the&e computations are shown in Figure 3 as a set of departure 
curves for the normal device. Curves are plotted for various values of Ru/Jt m 
(ratio of horizontal formation resistivity to mud resistivity) for four degrees of 
anisoiropy (X = l r 1.5,2, and.2.5). Each curve shows ihe variation of (ratio 
ot" apparent resistivity to mud resistivity) versus AM/d (ratio of normal spacing 
to hole diameter). 

The comments on an infinite, homogeneous, anisotropic medium are seen to 
describe the qualitative features of the curves quite accurately. For a given 
Rfi/R-,. increased anisotrcpy 'increased X) corresponds in general to an increased 
RJRm. (This explains why the 16-inch Normal reading at the level of shales may 
be higher than the Induction Log reading, as mentioned in the Introduction.) For 
spacing* large in comparison with Lhe hole diameter (AM/d large), the apparent 
resistivity j? a approaches Rh. 

It was possible to use this same approach to compute the effect of anisotropy 
when the formation has been invaded with mud filtrate. The formation close to 
the borehole will have a different resistiviiy from the uncontaminated part of the 
formation further back from the borehole. The formulas are then, of course, more 
complicated, particularly if the coefficient of anisotropy in the invaded zone is not 
assumed to he the same as in the uncontaminated zone. 

GENERAL CASE NXGLECTItfG TEE BOREHOLE AN"D INVASION 
If the influence of the borehole is neglected, the transformation to an isotropic 
medium can be made by setting 

/(p) - 1 and t(z) m X <2S) 

in equation (fin The transformation of equations (8) and (13) then involves only 
a stretching in the : -direction and no change in the radial distance or the mean 
resistivity. Each anisotropic bed *-ill be replaced by an isotropic bed having the 
same mean resistivity, but the thickness of the bed will be multiplied by the 
snisorropy coefficient X. Also, of course, after such a transformation, one muss 
lake for tbe new spacings between the electrodes of a sonde the differences in 
their coordinates, as given by fS>. 

The above transformation makes it possible to reduce the problem of andin? 
the potentials, or potential difierc.iccs, to be observed on the measuring elec- 
trodes ua the anisotropic case to that of ending the corresponding potentials in 
an equivalent isotropic case. To express the result* in terms of the apparent 
resistivity f? # recorded on a log, however, use is made of rhe linear relation 

R B = KV/I 
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FORUATIOS ASISOTSHOPY EFFECTS OS RESISTIVITY MEASVJt&lf&TS 
^d) ■=» V, at » *» «. 

where 



t = Vp 1 + s 1 and a is the borehole radius. 

Condition (a) speciaes the source as an isotropic current emitter located at 
the center of the borehole. Condition (b) simply fixes [he potential at infinity. 
Condition (c) is the statement of continuity of the normal component of current 
density across the cylinder representing the wall of the borehole. Condition 
(d) is a statement of the continuity of potential across the cylinder p—o and is 
equivalent to the statement thai there is no charge double layer on this surface. 

A formal solution ot this problem can be constructed in a fashion entirely 
analogous to the treatment given in the Appendix ot the Schlumherger Corpora- 
lion's Interpretation Bond-Book for Resistivity L^s. The notation adopted is that 
in the previous section, but now the origin is at the current electrode, and thus z is 
the spacing of the normal device. The potential in the borehole, expressed in terms 
of the modified Bessel functions Kn la, and I u is found to be 

Vm(p> t ) m -^[y + y f~A( X )h(z P /a) cos (xx/<) (22) 

where 

aW £ 

The "apparent" resistivity for tbc normal device is given by the formula 1 

R. - — V{Q t i). (23) 

The essential step in the computation of (23) is the evaluation of the Fourier cosine 
integral which appears in (23) when Vf.Q, t) U inserted from (22), namely: 

(2-4) 



p(\, A , i.) = 1 J%<*> cos (xz/a) dx. 



An effective technique lor the numerical evaluation of such an integral for 
large values of */<t and a slowly varying function 4<x) has been described by 
U N, G. Hlon (1928-1929). 

Filoa's suggestion amounts to filling A{x) with a set ot parabolic arcs- The 

v Tfc,. i MtUst a*f » the X. Of the normal Ocvkt. S« cquauon 1*o> «kj ww^tc^ 
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fORslATtOS AyiSOTROPY EFFECTS O.V RESISTIVITY UEASURZMESTS 

between And the measured* po cental difference, V. The proportionality con- ^/ 
stan: K is so chosen that R* is equal to the true resistivity for an intiniie, homo- 
geneous, isotropic medium. For the normal sonde, K is proportional to the spac- 
ing AM, whereas for the Lateral it is proportional to {A AT) ■ (A A" /{.VjV*>. Thus, 
if one computes directly the potential difference to be expected in any given 
case, by solving the equivalent isotropic problem above, this potential difference 
may be substituted for V in equation (2j) to obtain R+. 

Id some cases the equivalent Isotropic problem may have already been solved, 
as on existing interpretation charts tfrora r «i s t 0 r network measurements, ex- 
periments, or computations], giving the effect ol electrode spacing, bed thickness, 
etc. Such data will ordir.axily be presented in terms of an apparent resistivity 

JV - TO (2?) 

where the coefficient K' depends oo the spacing* used in ike equivalent isotropic 
problem. 

It is necessary, in using: such data, to know whether the K' of the equivalent 
problem is the same as the K in the original problem or whe:her a transforma- 
tion oi the K is involved. In the present case (in which the ^-co-ordinates, and 
hence the electrode spacing*, are transformed) there has indeed been a transfor- 
mation of the K. 

The apparent resistivity data of the equivalent case can be used to obtain the . 
potential difference V needed in (22) ii care is taken that at each depth the K* '^Jf 
must correspond to the transformed spaeings. Substituting this V into (22) r it is 
seen that the apparent resistivity (for the given anisotropic problem) is related 
to the apparent resistivity Jt,' (for the equivalent isotropic problem) through tbe 
equation 

r. - » r:. cm) 

Note that K is determined by the actual sonde spacings in the original problem, 
whereas K* is determined by the transformed spacing* in terms of. the ='-<o- 
ordinates. 

In particular, for the normal sonde, in the case being discussed of no borehole 
and no invasion, 




where is the average oi X between the current electrode A at and the 

measure electrode M at s*>z*. 

* For the aorn»l wnie tbe p<x«mial oiSewnce referred 10 ii ii*t between tbe electrode 5/ and * "--J 
rrfereact electrode H at the surface. See ta« SthJcurttrger Corporitioo's tnxrprclottom Band -Sink 
for Rtsvtm*y Lou, s>. 8-11. 
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For a lateral sonde having a short Af\Y span, the K' coefficient is very neariy 
proportional 10 {AOy/M A\ as measured in the equivalent isotropic medium, and, 
Therefore, 

1 X 

«(5o) - jg> -30} 



•where X* gi,zo) is the value of the anisotropy coefficient at Q, and X Li its average 
value, between A and 0. 

QUANTITATIVE STUDY Of THE TWO-MEDIA PROBLEM WITHOUT SOfcESOLE 

The problem Considered here belongs to the class treated in tie above section. It 
is assumed that two thick anisotropic beds have a common boundary, and i: is 
desired to evaluate the response oi the normal and lateral device* near this 
boundary. In the first bed the horizontal aad vertical resistivities are Ra % *nd 
/!*,, respectively . The m ean resistivity, = y/&* x • Rir) and the anisotropy co- 
efficient Xi( » VJZirJfcril are defined in agreement with equation (4). With 
similar dennitioos for the second bed, :be original problem can be replaced by 
an equivalent isotropic problem, in the manner indicated in the previous section 
\e-quaiion (25;). 

The current electrode A and a point of measurement M are assumed to be 
located on the same normal to the boundary between the two media. The upper 
hail of Figure 4 shows three different versions of the problem: (a) Electrodes A 
and Jf are in medium no. 1 ; ibi electrode .4 is in medium no. 1, and electrode M 
is in medium no. 2\ or (c) both electrodes are in medium no. 2- fin each case 
shown, A is above if. However, the same end result is found for the potential ii 
A and M are simply interchanged in position.) The equivalent isotropic cases are 
pictured in the lower half of Figure 4, 

U*ing L as the distance from electrode A to the boundary and : as the veit-ca! 
co-ordinate of Af with reference to A (positive downward), the potential in' the 
two beds can be found from the equivalent isotropic cases, by the method o: 
images. 

y = U!± [J— -l ~ ~( r 5 -Yj for A and Af in medium no. 1 ! 

4»\iL| = i Rx+R % \ . lL-s\ /J 

Y _ L . J** , 1 r SfiftS.?- 1 ; w 

4x R^R, i liL + X,(i - L) | : 
y = I2l.\JL+ Ry ~ Rt L ! -Y| for A and M in medium no. 2 ; 

Using iSece relailons, ihe apparent resistivities which would be fiven by 
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normal and lateral devices may be found.' For litis determination, the co- 
ordinal* ; is set equal to the spacing of the appropriate normal device. 

S pedal Cue A- Ri 

This is a. case of special interest because the equivalent isotropic problem cor. 
responds t.o a homogeneous medium. The Apparent resistivity profiles versus 
depth for such a. case are illustrated by the solid curves, In Fijjuxe 5 for the normal 
device, and in Figure 6 for the lateral. In each figure the upper anisotropic bed 
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Fie. 4. Two semi-infinite inbetroptc nedia, borehole effect acceded. 
CDrawn tor X,- 1.25. X»-2.) 



(region 1) has a horizon tat resistivity J?*, of 2, and the lower anisotropic bed 
(region 2) has Ra T equal to 1. (Note that the reference points for depths have been 
shifted to the midpoint of AM for ihe normal, and to point 0 for the lateral.) 

» The cormal apparent resistivity H computed Tram equation (2J> (<x equivslemiy equation <ZO)). 
In the present sad foJlo-winJ »eeiiona, the LauraJ appsreat resistivity is computed irwi tfee nili'ion 

!*t«r»J 

/ to 

whore > is set equaJ to the spaesog* AO. This rtladoo is s good appnndandoo wfcoj AfrY is negligible 
ia wb9»Hj«i with .40. 



S-t4 
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Also sbowc for comparison in Figures 5 and 6, in dashed line, are the apparent 
resistivity pronles for two isotropic beds whose resistivities R-* (upper bed) and 
J? j* (lower bed) are equal to the horizontal resistivities of tic corresponding 
anisotropic beds; i.e., and JV = \&* and Rf are not equivalent 

isotropic resistivities.) 
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Fio. 5. NontiAl wide, borcixoie effect a^kcied. 
Both normal and lateral curves for the anisotropic case show values oi R* 
which are constant and equal to ftAi), until the measure electrodes ctojj 

the boundary. Then, for the normal, while the A. electrode is in on* medium ar^ 
the measure electrode in the other, the value of the potential makes a smoo:!: 
transition from IRrJ*** to IJtiJtota. This transition from a» apparent resistiv- 
ity Rtr M to one of is completed in the distance z{-AM). which is one son.^ 
spacing- 

For the lateral, when the measure electrodes cross the bed boundary, i.-.e 

ipp^rent rewsuvitjr m»t*A & jump sv«W that the reading on the t*0 Side* 

the boundary are in the ratio X v % {equal in this case to RaJR» t ). This succc- 

J-13 
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jump is followed by a smooth transition t& tbc value R Hl { » which takes 

place in the distance of one sonde spacing AO. 

This example illustrates some of the novelty due to anisotropy. If we at- 
tempted to interpret the apparent resistivity curve in this case as though the 
media ^ere isotropic, we would have to assign resistivities and R Ht to match 




TK. & LaunJ icq it, borehole neglecled. 

the curve at a distance from the boundary. Then the behavior at the boundary 
would be misleading. (Compare dashed and solid curves in Figures 5 and 6.) 

Sfieeiei Case £ Ri*Rt 

The apparent resistivity it* will be the same in both media «cept ia the 
immediate neighborhood o£ the boundary. Let us suppose that X*>X|. Then, as 
indicated In figure 7, the reading of the normal device will increase as the 
boundary U approached from above, reaching a value 

\, + 
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when the measure electrode arrives at the boundary. Aher the measure elect rod? 
crosses the boundary, the apparent resisu'viiy depresses smoothly, becoming 
equal to 2Ai£«/(&i-r*ti when the current electrode reaches the boundary. The 
reading then increases to the value R* as the norma! device continues to move 
downward. 

The lateral device under the same ci:curr.»:zncci» will show a decreasing rtad- 
R 0 
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Pre. 7. Normal soode, borehole etJrti cwgteeteti- SpecoJ cue of R*. 
{Plotted forX,-l. )vi-2. ,8*«2,) 

iog. as shown in Figure S, arriving at a value 2Xi/2p/(Xi+X s ) when the measure 
electrodes are just above the boundary. When, the measure electrodes cross ihe 
boundary, the apparent resistivity jumps up by a factor of (\*Ai) 4 . ^ 
creases, reaching a. second minimum when (he current electrode passes the 
boundary. (The minima have the *amc value as the one for the normal device. 
Thereafter the apparent resistivity increases toward its £nal value of R*. 

For h/»th of the special cases illustrated in Figures 5 to 8 inclusive, it hjis be*:." 
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assumed that Xt>Xi, In a, case where Xt<A 5j the character of the features would 
be reversed. For example, in Figure S die- dip in the lateral just above the upper 
bed bound* ty would become a peak. The peak: just below the upper bed boundary 
would become a dip. The dip at one sonde spacing bciow the bed boundary wouid 
become a peak. 

The resistivity profile* of Figures 5 to S furthermore indicate that the response 
ft a - 

O t 2 * a 



A* --, 




?IC 8. Lateral sonde, boithtifc effect ftqdcCtcd- Sped*! case of J?*, X Mt fa. 
(plotted tor X t -2, X a mi.y 

of the lateral device is more strikingly aSected by the presence of anisotropy thin 
the response of the normal. 

General Case 

The most interesting features on the apparent resistivity curve air those ob- 
served when the measure electrodes are near the boundary and when the current 
electrode pr.ruz'z zic* the boundary. These value* are listed below m terms of the 
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"mean" resistivity 5, defined by 



R 2\R X Rj 



23 



measure electrode at boundary 



midpoint of AM it boundary 



NORMAL DEVICE 




current electrode at boundary 



(32) 



measure electrodes just above boundary 



measure electrodes just below boundary 



LATERAL DEVICE 



current electrode at boundary 



With these key values, a response curve may be sketched for any case of interest. 
If medium 2 is isotropic and we know Jt n , and fta from the readings above and 
below the boundary, then it is clear that readings corresponding to the lateral 
device response would allow the determination of Kv,(=*Xt.Ri) from the dis- 
continuity In &. at the boundary. 

The fact that the later*! device shows a discontinuity upon crossing a bound- 
ary in the ratio R^/Rf, is quite a general result. It is, of course, dependent ca 
the assumed s mattress of the Af\V-span between the measuring electrodes and cbe 
absence of borehole and invasion. 

THIN" tSOT=0t>TC $fl£© BZTWEEV ANISOTROPIC ADJACEVT FORMATIONS 'SAND 

oxr> aotnt dsd by shales) 
Using the technique o! reducing the anisotropic problem to an cquivalcss; 
isotropic problem, apparent resistivity profiles can be computed for the norma; 
and lateral devices passing through isotropic beds, such as sands, bounded by 
anisotropic adjacent formations, such as ahale*. From these results, interpreta- 
tion charts can be prepared. Only the lateral device will be considered here, be* 
causft it seems dear that this device will be influenced more significantly by 
anisotropy- 

That anisotropy in the adjacent formations may cause some ambiguity Lr. 
determining the resistivity of a thin bed is apparent from the fact that at eUfcer 

bou&d&fy the lateral reading in th« bad jusi inriA* the boundary, is relv.*-;'. 
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to the reading R*s. »f> the adjacent formation just ovixid* the same boundary, by 
the relation 

R.t Ra Rg 

where /?»• and }t„ are the vertical and horizontal resistivities In the adjacent forma- 
tions, X is the coefficient .of anisotropy in the 2djacent formations, and Rg is tbe 
resistivity ot the isotropic bed. At points wel] above or below* the bed the apparent 
resistivity is R+=Ru, If R a <>> i Rfr, equation (33) shows that jf? fl will decrease 
upon entering the bed. and this can occur even though Rb>Rh, Thus it is quite 
conceivable that a drop in £. wilt correlate with an increase la actual (horizontal; 
resistivity. From (33) this can be expected to happen /or 

l/,V < .3 < !. 

One might hope that this difficulty could be overcome by using other informa- 
tion from the resistivity curve, along with interpretation charts, to arrive at a 
better estimate oi Ra than the one deduced directly from the log. This possibility 
•**ill be considered. 

Special Com \Rjt = Rb 

Before considering the general case, however, let us take a look at a special 
problem that can be solved very simply by using the scaling technique of the 
sections which covered tbe general case (neglecting the borehole and invasion! 
and the quantitative study of the two-media problem without borehole. This 
special problem corresponds to the assumption that the mean resistivity of the 
anisotropic adjacent formation R t ( = y/'R v R H = \£ff)is equal to the bed resistivity 
R 9 . Tbe scaled problem then leads to a homogeneous medium of resistivity R*. 
and, therefore, it can be solved easily. This approach leads to an apparent 
resistivity protile for the lateral sonde as sketched En Figure 9. In this figure, s 
is taken equal to the lateral spacing .40. 

It is interesting to note that this curve could have been drawn directly from 
;2oj and (30), since In this case R& « Ra and the character of the curve is due 
completely to the variation of K'. Tbe curve obtained ir> this case is typical of 
cne result when. .S< 1, which is the case most likely to cause difficulty in log mtcr- 
preiation. II j5> 1 the discontinuities are executed in the opposite directions. 

General Cose *Rv?6R* 

In generalising to the case Rg^Rg, h i* important to decide what data from 
the Ra curve should b* used to determine Rg. It seems that for thin beda the best 
measurements would be: 

(1) a measurement at the upper or lower boundary of the ratio $ denned in 
(33); 

(2) a measurement of R* so the plauau when tbe sonde spans tbe bed- This 
value will be gWca tn* symbol & r - r 



7.10 
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T-.a. 9. LatermJ accdc, sotfbo^ ti«t argiected. Thin isotropic bed txtw«n ani»trop;t 
sioutcrr fomiiMnu tor special cav of Xm2.) 



(}) a measurement of R* in the adjacent formations; 
(4J a measurement of c, the bed thickness. 

By (33), RB/Rn^K&y vhere 3 is kno*'tt from measurement (I) above. Since 
in practice l<X<2.5 r it is known that p<Rg/Rfi <&.2S$. Far more precis 
knowledge of Rg, measurement (2) must be tnierprered using the faUo^'iftg 
theoretical expression, derivable by means of the theory of image*: 



where 



Re \2e/ \ 2e 2 J 

1 

a 1 — sceT^ 



a - the lateral spacing AO, 



i-ai 
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Using (33i and the definition of ct to eliminate X and « t in (34), it is possible to 
6nd aa expression for Rr/R/r which is a function of z/e, 0, and R a / Rh- Such a 
relation may be used to determine R 3 /Rh, since the other quantities are known 
from the four measurement* above. 

The nature oi the above so\ution is shown in the plots of R»/R rt vs R r /R H in ■ 
Figure 10a for z/c = l and in Figure 10b for r/e~5. Curves are plotted for sev- 
eral values of 3 (heavy curves). Alio shown are curve* of constant X (thin curves). 

Inspection of these curves shows that, particularly for the ;hin-bcd case 
= 5) for a given value of $, a small error in Rp/Rjt can cause a large error in 
Ra/R*. This means that for thin beds the method does not have much resolution 
for the determination of R B . For bed thicknesses approaching the lateral spacing 
the resolution is somewhat better (compare case of :/e=2). 

The above discussion cannot, of course, be taken to apply unmodified to the 
practical case where there is a borehole. U does point up the conclusion that there 
is no highly satisfactory way of correciing the lateral reading by itself to give 
the true resistivity of a thin isotropic bed between anisotropic adjacent iorrtW' 
lions. 



A method oE analysis has been described whereby a potential problem involv- 
ing anisotropic media can be transformed into an equivalent problem involving 
isotropic media only. The method is subject to the limitation of equation (6~> 
which requires ;hat the coefficient of anisotropy \ be expressible as a product of a 
function of the radial co-ordinate only by a function of the axial co-ordinate only. 

The potentials found by the solution of the equivalent isotropic problem, 
transformed back into the coordinates of the original problem, are sufficient to 
determine the apparent resistivities which would be observed by electrode* type 
logging devices In anisotropic formations. 

By means of the above transformations, one can piredtct, qualitatively, that 
a s.hor i -spacing measuring device in a borehole filled with isotropic mud will 
read higher than for an isotropic formation of resistivity R s (the horizontal 
anisotropic resistivity). This is confirmed to the case of the normal device by the 
accurately computed departure curves of Figure 3. 

In the cojc inhere the borehole if neglected some conclusions related to log inter* 
pre tat ion have been noted. 

1. Both normal and lateral devices measure Ru when the sonde is com- 
pletely within the bed, ana oot too close to a hed boundary. 



and 



Rb +■ Rs 




COXCIVHOV 
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(0) BED TVilCXNESS EOl^L CNE-HALF 



(b) BED THICKNESS EQUAL CNE-FlF" ru 
THE S»C3N<3 l*/e = 3» 



Fca. 10. L-itcril ton4c, borcbalc eSoct neglected, rcUUoa between Ra/Rtt » 
fuoexfca of d <"b«vy c^rv**) «ad X (ligM curve*). 



2, Seme anomalies occur when the electrode system is crossing a boundary 
between beds. In general these anomalies tend to be more striking for the lateral 
device. 

3. The apparent resistivity shown by the lateral device shows a discontinu- 
ity when the measure electrodes cross the bed boundary. The apparent resistiv- 
ities just above and just below the boundary are in proportion to the rcrt-cc/ 
Ttsiuivititi at these points (equations (32) and (33)). 

A. It b (lias pos&ibfc in sonc c**cs, lor the »snsc <?f tk« discontinuity in the 
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lateral curve to be in the opposite direction to the change in horizontal re- 
sistivity. As a result, a direct interpretation of the fog may sometimes be mis- 
leading particularly in the case of a thin isotropic bed between thick aniso- 
tropic shoulder formations. An attempt to interpret such a case from the read- 
ings of the later*! afonc indicated poor resolution in the determination of Aj, 
particularly for smaller bed thicknesses. 

In case there is a borehole, the above conclusions can apply approximately 
ii the aonde spacing is large in comparison with the borehole diameter. 
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Uet or Svw*M5 

Rm rvnnatUm resistivity weaaured by paasing current parallcJ to bedding plana (hofiiontal «t 
longitudinal ra&iiviry). 

Rr Formation resistivity measured by paasing correct in a direction perpendicular to bedding 
pla.ne* (verricjJ of trans versal t qatjwty), 

R Mean reactivity R*) : 

X Coefficient of anisotropy (X — %/R.r/Rn). 

p Radial coordlnaia of a cylindrical system of coordfn*te9- 

1 Axial cwdio*te a cylindrical system of ?5ordi nates, Ip Uie discussion, s is set equal 10 th« 

spacing of a normal or lateral reaislivity lodging device. \ ? 

$>' Radial coordinate in equivalent isotropic problem. 
%' Ajaa! coordinate in equivalent isotropic problem. 
R' Tbe fonnatioo resisrivity b the equivalent isotropic problem. 
Rm The resistivity of the isotropic mud colcoin filling the borehole 
A, Tbe apparent resisdvily indicated by 1 given resisuvixy logging device. 
R/ Tbe apparent resistivity »bfcb would be seen by the equivalent resistivity legging device la the 
tquivnk^t isotropic problem. 
R K reference resistivity for tbe two-media case, (R-7R*RJ(R,+RJ). wber« X, tube wean 

resistivity of the upper bed and K, is the nJcan tesisciviry oi the towrr bed. 
Ro Resistivity ot aa Uouopic bed. 

Ra Mean resistivity of anisotropic adjacent forma lions on either side of isotropic bed- 
R^t Apparent resistivity observed in a bed just within bed boundary. 

Apparent resistivity observed in anisotropic adjacent formation just outside bed boundary. 
jS Equal to R*s/R*f at either bed boundary In tbe no-borehole case. 
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